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We searched for the appearance of tt + tv~ pairs with invariant mass >648 MeV/c 2 in a neutral 
beam. Such an observation could signify the decay of a long-lived light neutral particle. We find no 
evidence for this decay. Our null result severely constrains the existence of an R° hadron, which is 
the lightest bound state of a gluon and a light gluino (gg), and thereby also the possibility of a light 
gluino. Depending on the photino mass, we exclude the R° in the mass and lifetime ranges of 1.2 - 
4.6 GeV/c 2 and 2 x 10" 10 -7 x 10" 4 s, respectively. 

PACS numbers: 13.85.Rm, 13.25.Es, 14.40.Aq, 14.80.Ly 

This Letter is motivated by recent discussions JlJ of the possible existence of long-lived hadrons that contain light 
gluinos. In some theories where the breaking of supersymmetry is communicated to ordinary particles by the exchange 
of very heavy states, gauginos have small tree-level masses and are light compared to the squarks. The gluino(g) and 
photino(7) masses are expected to be < 1.0 GeV/c 2 . The gluons (g), gluinos, and quarks can form bound states, 
the lightest of which is a spin-1/2 gg combination called the R°. The approximate degeneracy of R° with the ++ 
gluonia [jlj suggests that the R° mass could be of the order of 1.3-2.2 GeV/c 2 . The stable photino in this theory is 
a cold dark matter candidate ||. Estimates from particle physics |ij and cosmology Q place the R° lifetime in the 
10 _10 -10 -5 s range. The R° decays into a 7 and hadrons. Because of the approximate C invariance of supersymmetric 
QCD, the R° decay into is expected to be the dominant decay mode 0; depending on the extent of C violation, 
the R° may also decay into 7r°7 or 777. The existence of a light gluino can have an impact on the running of the 
strong coupling constant a s . A phenomenological analysis of the perturbative running of a s ||, in conjunction with 
the multijet analysis at LEP [1 has claimed an indirect exclusion of the light gluino scenario. However, a debate over 
the extent of this exclusion |5| underlines the necessity of a direct search for hadrons containing light gluinos. 

The i?°'s can be produced in pN collisions by processes such as quark-antiquark annihilation, gluon fusion, etc. 
Since squarks need not be involved, the R° production is in the realm of conventional QCD. Dawson, Eichten, and 
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Quigg |g[ have calculated cross sections for gluino production in tree approximation. Their calculations suggest an 
order-of-magnitude cross section estimate of ~ 10 /ib per nucleon for the production of an R° with 2 GeV/c 2 mass 
in 800 GeV/c pN collisions. This cross section corresponds to ~ 10 -3 R° to Kl flux ratio in our experiment, as 
explained later. Assuming a 10 -8 s lifetime and a 100% branching ratio for the decay of this R° into jp, p — > 7r + 7r~, 
the KTeV experiment Q has a per spill sensitivity of approximately 10 -4 in terms of the R° to Kl flux ratio with 
3.5 x 10 12 800 GeV/c protons delivered in a nominal spill. The results presented here are from the delivery of 1.9 x 10 15 
protons. 

In the KTeV experiment (figure |l|), protons are incident on a 30 cm long beryllium oxide target at a vertical angle 
of 4.8 mr with respect to the neutral beam channel. The interaction products are filtered through a 50.8 cm beryllium 
absorber and a 7.6 cm lead absorber. Two neutral beams, each 0.25 /istr in solid angle, emerge after collimation and 
sweeping. One of the beams passes through an active regenerator, but the decays from this beam are not used in this 
analysis. The beam transport and decays took place in an evacuated region with a vacuum of 0.5-1.0 x 10~ 4 torr. 
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FIG. 1. KTeV detector configuration for this measurement. Note the highly compressed scale along the beam direction. 



The most crucial detector element for this analysis, the charged spectrometer, consists of four planar drift chambers, 
two on either side of an analyzing magnet. Each chamber measures positions in two orthogonal views. Each view 
consists of two planes of wires, with cells arranged in a hexagonal geometry. Each chamber has approximately 100 /im 
single-hit position resolution per plane. The spectrometer magnet has a 2 m by 1.7 m fiducial region over which a 
transverse momentum impulse of 411 MeV/c is imparted to the charged particles. A set of helium bags integrated 
into the spectrometer system minimizes multiple scattering. The invariant mass resolution for the decay Kl — ■> 7r + 7r _ 
is better than 2 MeV/c 2 . 

A pure Csl electromagnetic calorimeter of dimension 1.9 m by 1.9 m is used to reconstruct photon and electron 
energies to better than 1% precision. The calorimeter is used to match the orthogonal track views and to reject 
background from K e 3 decays. A set of 12 photon vetos provides hermetic photon coverage up to angles of 100 mr. 
A counter bank (muon veto) located at the downstream end of the detector is used to reject decays. The event 
trigger is initiated by signals from two scintillator hodoscopes located downstream of the spectrometer. The primary 
trigger requires two hits in these counters consistent with two oppositely charged tracks, at least one hit in each of the 
two upstream drift chambers, and lack of hits in the muon veto bank. Next, a set of fast trigger processors requires 
two straight tracks in the non-bend view from two hits in each drift chamber. 

A combination of online and offline cuts in the 7r + 7r~ analysis required two oppositely charged tracks matched to 
clusters in the calorimeter to within 4.6 cm. Individual track momenta were required to be more than 8 GeV/c, and 
the scalar sum of two track momenta was required to be between 30 GeV/c and 160 GeV/c. A longitudinal vertex 
position between 126 m and 155 m downstream of the target defined the fiducial region. Electrons were rejected by 
requiring the energy deposited in the calorimeter to be less than 80% of the particle momentum, as measured in the 
spectrometer. The signals in various veto devices were required to be no more than those due to accidental activity. 
The two tracks were required to project within the active fiducial area of the muon banks. After making all other 
cuts, the momentum ratio of the two tracks was required to be between 0.2 and 5, which curtailed the high-side tail 
of the 7r + 7r~ mass distribution by moving its end point from ~640 MeV/c 2 to ^600 MeV/c 2 . Figure || shows the 
m w + w - distribution for the events surviving all cuts. Note the peak corresponding to the Kl — > 7r + 7r~ candidates 
and the rapidly falling background to the right of the kaon peak. The figure also shows the m 7r + 7r - distribution for 
an R° with mass m R o=1.75 GeV/c 2 and photino mass m^=0.8 GeV/c 2 . There are no R° candidates in the signal 
window between 648 MeV/c 2 (i.e. 150 MeV/c 2 above the kaon mass) and 1.0 GeV/c 2 . For the simulated R° shown, 
91% of the decays are contained in the signal window. 
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FIG. 2. 771,,.+^- distribution for the data (solid) and R° signal Monte Carlo (dashed, arbitrary scale). The peak at 500 MeV/c 2 
corresponds to Kl — > tt + tt~ decays. The sharp cutoff at 0.95 GeV/c 2 (=m R o — rrij) is due to the kinematic limit for the R° 
shown. 



The absence of signal in the data can be expressed in terms of upper limits on the R° flux. We define the R° /Kl 
flux ratio to be the ratio of the number of R°'s to the number of K^s exiting the beam absorbers, calculated with 



the assumptions that the R 



IP, P 



"branching ratio is 100% and that the photino does not interact 



significantly in the detector material. For the R° spectrum shape, we use the invariant cross section for A production 
in pBe interactions ||. We make this choice because the R° mass is in the proximity of the A-S°-Z?° mass range. 
For our experiment, the function (1 — ccf) exp(— bp\) with a = 1.0 and b = 2.3(GeV/c) -2 is a good approximation 
for the cross section shapes given in |^] for both A and S°. The limits are not very sensitive to the spectrum shape; 
if we use the values a = 6.1 and b = 1.08(GeV/c)~ 2 which are applicable to D° production ]l0|] , the limits change by 
-5-10% for R° lifetime of 10" 8 s. 



The Kl flux was determined using a total of 116,552 Kl 



candidates with two-body transverse momentum 



squared (pj) less than 250 MeV /c and m 7r + 7r - within 10 MeV/c of the Kl mass. The detector acceptance for the 
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decays was calculated using a Monte Carlo simulation of the beam and the detector. We calculate a 
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FIG. 3. Geometrical acceptance of the detector in the fiducial region for various R° masses as a function of the R° momentum. 
m R a and are in GeV/c 2 . (R° lifetime = 10 _8 s) 
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total of 1.33 x 10 10 Kl's of all energies exiting the absorbers. To calculate the detector acceptance for R 0, s, the 
R — > 7P, p — > 7r + 7r _ decays were simulated for various values of m^jo and assuming isotropic angular distributions 
of the decay products in the center of mass frame. Figure || shows the geometric acceptance of the detector in the 
fiducial region for various R° masses as a function of the R° momentum. Figure ^ shows the 90% confidence level 
upper limits on the R°/Kl flux ratio as a function of m^o with a 10~ 8 s R° lifetime for two values of the mass 
ratio r = m^o/m^. The sharp loss of sensitivity for small m^o for a given r is because an R° with mass less than 
0.648r/(r — 1) GeV/c 2 can not produce a 7 with mass m R o jr together with a 7r + 7r~ pair having invariant mass greater 
than 648 MeV/c 2 . 
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FIG. 4. Upper limits with 90% confidence level on the R°/Kl flux ratio as a function of R° mass for two different values of 
the i?°-photino mass ratio r. (R° lifetime = 10~ 8 s) 

Figure |B| shows the upper limits on R°/K~l flux ratio as a function of R° lifetime with different masses for the same 
r(=2.2). Figure ^ shows the 90% confidence level upper limit contours for various R° masses and lifetimes. These 
upper limits can also be expressed in terms of limits on the invariant cross section Ed 3 a /dp 3 times the branching 
ratio using the conversion factor of 2.1 x 10~ 36 cm 2 /(GeV 2 /c 3 ) at xf — 0.2 per 1 x 10~ 7 R°/K~l flux ratio . 
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FIG. 5. Upper limits with 90% confidence level on the R° /Kl flux ratio as a function of R° lifetime. m R o and are in 
GeV/c 2 .(r=2.2) 
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Our limits are a significant improvement over the previous (indirect) search by Bernstein et al. |Tl[. Note that this 
conversion factor and the upper limit on Ed 3 a /dp 3 reported by Bernstein et al. do not take into account the R° 
absorption in the target and absorbers. The R°N cross section is expected to be in the range 1/10 to 1 times 

the NN cross section, which places the R° absorption factor in the 1.3 to 10.2 range. 




FIG. 6. Upper limits with 90% confidence level on the R°/Kl flux ratio for various R° masses and lifetimes. (r =2.2) 

The upper limits on the RP /Kl flux ratio constrain the light gluino scenario assuming a specific model for the R° 
production. The perturbative QCD calculations suggest a(R°) ~ 860e _2 - 2m K° ^b in 800 GeV/c pN collisions, 
where m R o is in the units of GeV/c 2 . This estimate is approximately consistent with the heavy flavor production 
cross sections. For our experiment, this <j(R°) implies that the R? /Kl flux ratio for an R° of mass mpo is expected 
to be 9.2 x 10 -2 e -2 ' 2m - R ° , assuming a factor of 10.2 for the R° absorption. The mass-lifetime region for which this 
R° I Kl flux ratio expectation exceeds the measured upper limits is taken to be ruled out. For example, the expected 
R° /Kl flux ratio for an R° with 2 GeV/c 2 mass is 1.1 x 10~ 3 . Since this expectation exceeds our upper limits in the 
lifetime range of 3.4 x 10 -10 - 1.3 x 10~ 4 s, we rule out R°'s with 2 GeV/c 2 mass in this lifetime range with 90% 
confidence level. Figure R shows the excluded regions obtained in this fashion for two different values of the mass 



p, 5 
o 




■| r i i i i i 




FIG. 7. R° mass-lifetime regions excluded by this analysis at 90% confidence level for two different values of the i?°-photino 
mass ratio r. 
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ratio r. Apart from the low mass cutoff for m^o < 0.648r/(r — f) GeV/c 2 , the contour shapes in figure fj] are fairly 
insensitive to the value of mass ratio r. The contours for the cosmologically interesting region correspond to r < 1.8 [0. 

In summary, our search places stringent upper limits on the R^/Kl flux ratio . For example, the 90% confidence 
level upper limit on the i?° /Kl flux ratio for an R° with mass 1.75 GeV/c 2 , lifetime 6 x 10 -9 s, and mass ratio r~2.2 
is 1.2 x 10~ 7 . As figure ^ shows, these stringent limits rule out the production of R° particle in a wide range of masses 
(~1.2 - 4.6 GeV/c 2 ) and lifetimes (~3x 10~ 10 - 7 x 1CP 4 s). A search |l3| for light supersymmetric baryons (uudg 
and udsg) in the mass range 1.7-2.5 GeV/c 2 has also reported null results. 
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